On the environmental stability of quantum chaotic ratchets 
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The transitory and stationary behavior of a quantum chaotic ratchet consisting of a biharmonic 
potential under the effect of different drivings in contact with a thermal environment is studied. 
For weak forcing and finite h, we identify a strong dependence of the current on the structure of 
the chaotic region. Moreover, we have determined the robustness of the current against thermal 
fluctuations in the very weak coupling regime. In the case of strong forcing, the current is determined 
by the shape of a chaotic attractor. In both cases the temperature quickly stabilizes the ratchet, 
but in the latter it also destroys the asymmetry responsible for the current generation. Finally, 
applications to isomerization reactions are discussed. 
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I. INTRODUCTION 

Directed transport, understood as transport phenom- 
ena in periodic systems out of equilibrium, has attracted 
much attention in recent years [1, 2]. As a consequence 
the field has developed into a well established area of 
statistical physics, which involves many interdisciplinary 
aspects. In this respect, many classical questions and 
very recently many quantum issues have been answered. 
The breaking of all spatiotcmporal symmetries leading 
to momentum inversion has been found to be the gen- 
eral mechanism to engineer ratchet systems [3]. For the 
Hamiltonian case, an efficient sum rule explaining the 
values of the resulting net current has been devised [4] . 
Quantum effects were considered to analyze the first so- 
called quantum ratchets [5], while recently purely quan- 
tum ratchets have been found to exist [6] . Floquet theory 
has provided with a general explanation for the appear- 
ance of a quantum current in periodic systems [7] . 

Ratchet systems are interesting in a very broad range 
of situations, such as in applications to molecular motors 
in biology [8] or to nanotechnology [9] . Cold atoms in op- 
tical lattices are one of the main examples of successful 
implementations and theoretical developments [10, 11]. 
Also, Bose-Einstein condensates have been transported 
(for particular initial conditions) by using purely quan- 
tum ratchet accelerators [12]. In this case, the current 
has no classical counterpart [13], and the energy grows 
ballistically [14, 15]. A possible application field that has 
not been much explored yet is represented by molecular 
processes, like for example isomerization. This particular 
type of chemical reactions have a tremendous relevance 
in important biological processes, such as human vision 
[16] or proton transfer [17], and their control has also 
been considered in the literature [18]. Actually, some of 



us have very recently proposed a novel method to per- 
form this control in a realistic model for LiNC^LiCN 
[19], which is an adequate prototypical example for this 
kind of processes. Notice that in these studies the main 
ingredients used in the present paper, namely dissipa- 
tivc dynamics associated with thermal noise and external 
perturbations [20], are included. However, they were not 
formulated and studied from the directed transport per- 
spective. Here we focus on the behavior of a dissipativc 
ratchet in regimes that can also be of interest for this 
kind of experiments. 

In this paper we study the influence of the environment 
on a quantum chaotic ratchet. Our model consists of a 
mass particle in a biharmonic potential subjected to dif- 
ferent periodic driving forces. In the dissipative case, the 
environment can be directly responsible for the transport 
generation. Let us remark that the results for this regime 
are still scarce since the calculations are difficult to carry 
out, and then many questions remain still unsolved [21- 
23]. New results seem to provide some solutions to these 
problems [24, 25]. Among them, determining the stabil- 
ity of the current is of great relevance and has recently 
raised high interest [26] . We show that for weak forcing 
and finite h values the current strongly depends on the 
structure of the chaotic region. Thermalization brings 
stability, without completely washing out the asymmet- 
ric structures that give rise to it. In the strong forcing 
scenario, the current is explained by the asymmetry of 
a chaotic attractor. A stable classical and quantum cur- 
rent is achieved at short times, and the effect of moderate 
temperature consists of making this times even shorter, 
but at the price of diminishing the current value. 

The organization of this paper is as follows. In Section 
II we present our model for the system and the envi- 
ronment, also describing the methods used to investigate 
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the current behavior. In Section III we show the results, 
and the roles of the coupling strength and the temper- 
ature are analyzed in detail. Finally, in Section IV we 
summarize our conclusions. 



II. MODEL AND METHODS 

Our system consists of a particle moving in a time de- 
pendent potential given by 

V(x,t) = 1 - cos(x) - A cos(2x + 4> a ) + (1) 
Fs'm(x) [cos(t) + Bcos(2t + &)], 

where F is the strength of the time periodic forcing, A 
and 4>a are parameters that allow to introduce a spatial 
asymmetry, and B and <pb have an analogous effect in 
the time domain. Throughout this paper we will take 
A = 1/2 and 4>b = tt/2. The effects of the environment 
are introduced by means of a velocity dependent damping 
and thermal fluctuations. This leads to integrate the 
equation 

mx = -rx-V'(x,t)+£. (2) 

In the above expressions x is the spatial coordinate of 
the particle, m its mass, and T, with < T < 1, is the 
dissipation parameter. The thermal noise £ is related 
to T, according to < £ 2 >= 2(1 — r)fcsT, where ks is 
the Boltzmann constant and T is the temperature, thus 
making the formulation consistent with the fluctuation- 
dissipation relationship. In the following we take m = 1 
and fee = 1- 

At the quantum level, we perform the evolution of the 
density matrix of the system, p, by means of a modi- 
fied split operator method [27]. We use a composition 
of unitary steps given by the kinetic and potential terms 
of the Hamiltonian (representing the system dynamics), 
and other purely dissipative steps. The latter come as 
the result of incorporating dissipation and thermaliza- 
tion to the quantum particle by coupling it to a bath 
of noninteracting oscillators in thermal equilibrium at a 
temperature T. The degrees of freedom of the bath are 
eliminated by means of the usual weak coupling, Markov 
and rotating wave approximations [28]. As a result, we 
arrive at a Lindblad equation for the density matrix of 
the system, that can be written as a completely positive 
map ~D a (dt) in the operator-sum (or Kraus) representa- 
tion 

2 

p(t + dt) = D (eiT) (tft) (p(t)) = ^P(*)^ f . ( 3 ) 

ft=0 



where 

i 2 

K = 1 - - £ K?K± 

N-X 

Kf = \A dt (1 + fi(fc)) k \p± kT i)(p±k\ 

k=l 
JV-1 

K± = J2 Vzdtn(k) k \p± k )( P±kTl \ (4) 

k=l 

arc the infinitesimal Kraus operators satisfying: 
Y,^ K P K ti = 1 to first order in dt [29], and p is 
the momentum conjugated to the x coordinate. Note 
that the superscript ± defines two different operators 
(standing for the positive and negative values of the p 
spectrum), and do not apply to the Kq operator. In 
these equations e is a system-bath coupling parameter 
that can be directly associated to the classical velocity 
dependent damping T (at T = gives the contraction 
rate of the phase space). The population densities of 
the bath are given by n = (exp(Ai5fc/(fcsT)) — l) -1 , 
where we have taken E k = p|/2, and AE k is the energy 
difference between the neighboring levels connected by 
the operators. In this way, we extend our method [23], 
originally developed for maps, to general fluxes. 

III. RESULTS 

In directed transport studies the main quantity that 
characterizes the system is the current J{t) = (pt), where 
() means the average taken with respect to the initial con- 
ditions and time, up to a given instant t. In the quantum 
case, we consider J(t) = (Tr(pp)), where the same kind 
of average is taken. We study its transitory and station- 
ary behavior. Given that phase space distributions are 
also of great interest, in the following we will also focus 
on their analysis. 

We have considered two cases, weak and strong forc- 
ing. In the first case (F ~ 0.02 — 0.05) we have bro- 
ken all spatiotemporal symmetries that forbid a net cur- 
rent (namely, a generalized parity and time reversal) by 
means of the potential built-in asymmetry (we have taken 
B = 0.2 and (j> a =0). In this situation the phase space 
of the closed system is mixed and the current depends 
on the structure of the chaotic sea, which has to be 
asymmetric. Moreover, we have tested the robustness 
of this mechanism against dissipation and thermal fluc- 
tuations. In the second case (F = 2.5) we have chosen 
a simpler, harmonic forcing {B = 0), and an asymmetric 
spatial potential (4> a = 7r/2) in order to investigate the 
interplay among forcing, dissipation, and thermal fluc- 
tuations. Now, the current comes from the asymmetry 
of a strange attractor. In fact, dissipation induces this 
asymmetry, which is responsible for the directed trans- 
port [21]. But this same dissipation mechanism contracts 
phase space and makes the higher energies inaccessible 
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for the system. So higher dissipation and asymmetry 
does not translate necessarily into higher values of J. 
Additionally, thermal noise compensates for the energy 
loss caused by dissipation. However, the associated diffu- 
sion tends to homogenize the strange attractor, reducing 
its asymmetry, and these two effects compete with each 
other [23]. 



A. Weak forcing 

The biharmonic spatial potential we have considered 
has two minima. For weak time-periodic forcing there is 
a chaotic region that connects two corresponding main 
regular islands which are related to these minima. In 
fact, a chaotic region fills the portion of phase space 
among many regular islands, in the shape of a network 
of branches, with almost no bulk portion. As a conse- 
quence, there is a high sensitivity to variations in the 
intensity of the driving field, since this introduces rele- 
vant modifications in this kind of ballistic network. This 
phenomenon can be clearly seen in Fig. 1, which displays 
the phase space for F = 0.02 (upper row) and F = 0.05 
(lower row) at time t = 50. In all cases the coupling with 
the environment is very weak, T = 10~ 4 . In addition 
we take h = 0.041, which is far from the semiclassical 
limit. The classical and quantum distributions are ob- 
tained by evolving a set of analogous initial conditions 
inside the chaotic region with initial p = 0. The contour 
lines indicate different level surfaces of the classical dis- 
tributions while the density plots show the quantum ones 
(where darker means higher values). The imbalance be- 



which is responsible for the current, is significantly al- 
tered by varying the parameter F. For F = 0.02 and 
T = (see Fig. 1, upper row, leftmost panel) there is 
a branch of the chaotic region that develops for p < 
and that is not present for p > 0. For F = 0.05 and 
T = (sec Fig. 1, lower row, leftmost panel) the distri- 
bution looks more symmetric. This behavior persists at 
low temperatures (T = 0.01) as can be seen in the mid- 
dle column of Fig. 1, for both values of F. Finally, for a 
higher temperature (T = 0.1) the situation changes, and 
the distributions look more alike, although some features 
of the phase space structure seem to survive the thermal 
fluctuations. 

In order to assert if this is the case, we show the cur- 
rent values as a function of time (in units of the period 
of the forcing) in Figs. 2 (a) (for F = 0.02) and (b) (for 
F = 0.05). We notice that at T = (blue dot-dashed 
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Figure 1: Classical (contour lines) and quantum (gray scale 
density plot) phase space distributions for t = 50 (in units of 
the period of the forcing). In the upper row F = 0.02 and in 
the lower one F — 0.05. We display the cases (from left to 
right) for T = 0, T = 0.01, and T = 0.1. In all cases Y = 10~ 4 
and h = 0.041. 
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Figure 2: (Color online) Quantum current as a function of 
time t (in units of the period of the forcing) for the cases 
shown in Fig. 1. In panel (a) F = 0.02 and in (b) F = 0.05. 
(Blue) dot-dashed lines correspond to T = 0, (orange) dashed 
lines to T = 0.01, and (black) solid lines to T = 0.1. In all 
cases T = 10~ 4 and h = 0.041. 



tween the positive and negative p regions of phase space, 



lines) there are strong fluctuations that can be associ- 
ated with the very weak coupling with the environment 



4 



and the finite H value. At T = 0.01 (orange dashed lines) 
the same behavior is present, but at T = 0.1 (black solid 
lines) the situation changes and the current stabilizes at 
very short times (t ~ 15), so that a value for the cur- 
rent J can be defined. Moreover, this is valid for both 
values of F. At this point it is worth mentioning that 
for F = 0.02 the current is greater than for F = 0.05, so 
we can conclude that the quantum distributions preserve 
the main features of the original structure of phase space 
at T = 0. This makes the current generation mechanism 
robust against environmental perturbations in the weak 
coupling regime. Moreover, temperature also helps to 
stabilize the ratchet. 

To conclude this subsection, we would like to men- 
tion that our choice of parameters is suitable for model- 
ing many isomerization reactions induced by laser fields 
[16, 18, 20]. The interactions with the solvent can ac- 
count for the interaction with a thermal environment of 
the kind that we have analyzed. In this respect, it is in- 
teresting to note that the main features for current gen- 
eration (isomerization rate in this case) survive a weak 
coupling with the environment, and that this in turn can 
be beneficial. This suggests that this mechanism is ap- 
plicable to actual experimental situations. 



B. Strong forcing 

The same system considered above can behave in a 
quite different way if the strength of the forcing is in- 
creased. In fact, when this happens the regular islands 
structure present in the previous case is completely lost, 
and then the current arises as a consequence of the asym- 
metry of a chaotic attractor. In this case, it is easier 
to define an asymptotic current given the fact that the 
attractor is usually formed in a very short time. To il- 
lustrate this effect, we present in Fig. 3 the classical and 
quantum phase space distributions for F = 2.5, T = 0.05, 
and H = 0.041 at t = 50. In it, we show from left to right 
and top to bottom the cases corresponding to: T = 0, 
T = 0.001, T = 0.01 and T = 0.1, respectively. As 
in the previous case, the influence of the temperature 
is noticeable only for its higher value, namely T = 0.1; 
while for the rest of the cases the quantum distribution 
remains almost unchanged. Thermal fluctuations need 
to be greater than the quantum ones to manifest them- 
selves. This is not the case for the classical attractor, 
which gradually loses all its finer details. For the set of 
parameters that we have chosen in this case, the time 
reversal invariance is only broken by dissipation. More- 
over, dissipation is not only responsible for the current 
generation, but also for the quick stabilization of the cur- 
rent. This can be seen in Fig. 4, where the classical (thin 
lines) and quantum J as a function of time t are shown. 
As can be seen, for T = and T = 0.01 the attrac- 
tor needs ~ 10 periods of the forcing to set in; at this 
point the current is very well defined, both classical and 
quantum mechanically Actually, the correspondence be- 
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Figure 3: Classical (contour lines) and quantum (gray scale 
density plot) phase space distributions for t = 50 (in units 
of the period of the forcing). From left to right and top to 
bottom we show the results corresponding to: T — 0, T = 
0.001, T = 0.01 and T = 0.1, respectively. In all cases F = 
2.5, T = 0.05 and ft = 0.041. 




Figure 4: (Color online) Classical (thin lines) and quantum 
current as a function of time t (in units of the period of the 
forcing) for some of the cases shown in Fig. 3. (Blue) dot- 
dashed lines correspond to T = 0, (orange) dashed lines to 
T — 0.01, and (black) solid lines to T = 0.1. In all cases 
F = 2.5, T = 0.05 and ft = 0.041. In the inset we show 
enlarged the stabilization region for the sake of clarity. 
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tween the quantum and classical results is remarkable. 
For T = 0.1 the setting in of the attractor is even faster 
(just ~ 6 periods are needed), but in this case the effect 
of temperature on the quantum current is rather drastic, 
i.e. J goes to zero. We have further verified this behavior 
by using a highly efficient diagonalization method [30]. 
This allowed us to look into the details of the superop- 
crator spectrum, for which the spectral gap is around 
0.4 (this gap is the distance between the modulus of the 
leading eigenvalue - equal to one - and the modulus cor- 
responding to the next eigenvalue, in descending order 
of this quantity) . We have also obtained the associated 
equilibrium eigenvector (corresponding to the eigenvalue 
one), which resulted indistinguishable from the chaotic 
attractor found by means of the time evolution. 

In order to show the ability of our method to reach the 
semiclassical limit we consider the classical and quantum 
chaotic attractors corresponding to h ~ 0.0015 and T = 
in Fig. 5. Let us mention here that the study of the 
evolution of different quantities as a function of time is 
an interesting sub-product of our approach. Thanks to 
this, we have found that despite the greater accuracy 
in the details of the phase space distributions attained 
in this case, the values of J as a function of t remain 
essentially the same as for h = 0.041. 




2ti X 

Figure 5: (Color online) Classical (left panel) and quantum 
(right panel) phase space distributions at time t — 50 cor- 
responding to T = 0. In all cases F — 2.5, T = 0.05 and 
H~ 0.0015. 



IV. CONCLUSIONS 

In this paper, we have studied a ratchet system con- 
sisting of a particle moving in a biharmonic potential 
subjected to a time periodic forcing and in contact with 
a heat bath at finite temperature. The rich structure of 
the phase space for different values of the involved pa- 
rameters has been analyzed both classical and quantum 
mechanically. 

We have found that for weak forcing and very weak 
coupling with the environment, the mechanism leading 
to current generation survives. In fact, the quantum dis- 
tributions keep "memory" of the original structure of the 
chaotic region at T = 0. Moreover the effect of the tem- 
perature is beneficial, rather than destructive, and the 
current stabilizes after short times. This makes us sug- 
gest a possible experimental verification of these results 
in chemical isomerization reactions, to which the dynam- 
ics we have analyzed can be directly applicable. In the 
case of the strong forcing regime, the effect of tempera- 
ture also shorten the times for which the stabilization of 
the current occurs. But in this case, the current vanishes 
for higher values of T, due to the blurring of the chaotic 
attractor details (and asymmetry). 

Finally, it is worth mentioning the high precision and 
simplicity of our integration method. By developing a 
modified split operator scheme, we have been able to 
reach very low values of %. This, combined with a 
highly efficient diagonalization method will allow us to 
implement in the near future an alternative approach 
for studying time dependent dissipative quantum systems 
[30]. 
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